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ABSTRACT 

We present a comparison between the observed properties of damped Lya sys- 
tems (DLAs) and the predictions of simple models for the evolution of present day 
disk galaxies, including both low and high surface brightness galaxies. We focus in par- 
ticular on the number density, column density distribution and gas density of DLAs, 
which have now been measured in relatively large samples of absorbers. From the com- 
parison we estimate the contribution of present day disk galaxies to the population 
of damped Lya systems, and how it varies with redshift. Based on the differences be- 
tween the models and the observations, we also speculate on the nature of the fraction 
of DLAs which apparently do not arise in disk galaxies. 

Key words: galaxies: evolution, galaxies: formation, galaxies: spiral, quasars: ab- 
sorption lines 



1 INTRODUCTION 

The true nature of the galaxies responsible for high red- 
shift damped Lya systems (DLAs, absorbers seen in quasar 
spectra with H I column densities N(U I) > 2 x 10 20 cm" 2 ) 
is largely unconstrained. At present, there are two main 
competing scenarios for their origin. One school of thought 
sees DLAs as the (large) progenitors of massive spiral disks 
(Wolfe et al. 1986; Lanzetta et al. 1991). The gas disks would 
have formed at z > 5 through monolithic collapse, and this 
gas is converted to stars over a Hubble time. In support of 
this picture Prochaska & Wolfe (1998) argued that the kine- 
matics of the metal absorption lines in DLAs are best ex- 
plained if they are formed in thick, large, and rapidly rotat- 
ing galactic disks, with circular velocities Vo ~ 200 km s . 
However, this large disk hypothesis runs counter to currently 
popular models of hierarchical structure formation in which 
present day galaxies are assembled from virialized sub-units 
over a protracted time interval [z ~ 1 — 5). Haehnelt et al. 
(1998) argued that hydrodynamic N-body simulations are 
able to reproduce the velocity structure of the absorption 
lines with infalling sub-galactic clumps in collapsing dark 
matter haloes with small virial velocities (V ~ 100 km s _1 ). 

More generally, the kinematics and metallicites of DLAs 
have been interpreted as evidence that they arise in spiral 
galaxies (Fritze-V.Alvensleben et al. 2001; Hou, Boissier & 
Prantzos 2001), low surface brightness galaxies (Jimenez, 
Bowen & Matteucci 1999), dwarf galaxies (Matteucci, Mo- 
laro & Vladilo 1997), the progenitors of globular clusters 



(Burgarella, Kissler-Patig, & Buat 2001), the building blocks 
of current galaxies (Tissera et al. 2001), galaxies undergoing 
tidal stripping and mergers (Mailer et al. 2001), and outflows 
from dwarf galaxies (Schaye 2001a). 

Observationally, imaging studies of the fields of QSOs 
with damped systems have shown conclusively that the ab- 
sorbers are a very 'mixed bag', which includes galaxies of 
different luminosities and surface brightnesses, down to ob- 
jects with apparently no associated stellar populations which 
remain undetected even in very deep images (Steidel et al. 
1994, 1995; Le Brun et al. 1997; Lanzetta et al. 1997; Fynbo 
et al. 1999; Pettini et al. 2000; Turnshek et al. 2001; Bowen, 
Tripp, & Jenkins 2001; Kulkarni et al. 2000, 2001; Colbert 
& Malkan 2002). While it seems clear that selection based 
on H I absorption cross-section picks out a variety of galax- 
ies, it is of interest to establish if one particular class of ob- 
jects dominates and, if so, whether the dominant population 
changes with cosmic epoch. 

In this paper we assess the contribution of the progeni- 
tors of today's disk galaxies to the DLA population at differ- 
ent redshifts by comparing the most up to date observational 
determinations of several properties of DLAs with models of 
the chemical and spectrophotometric evolution of disk galax- 
ies, including low surface brightness (LSB) galaxies. Specif- 
ically, we have compiled recent data on the number density 
per unit redshift, the column density distribution and its 
integral, which gives the total mass of H I traced by DLA, 
for a large sample of DLAs. These data are presented in §2, 
while the models of disk galaxies are discussed in §3. In §4 
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we compare our predictions with the statistical properties 
of DLAs, as well as with available imaging observations of 
the galaxies identified as DLA absorbers. We further specu- 
late on what the differences we find at high redshift between 
models and observations may be telling us about the earliest 
population of DLAs. We summarise and discuss our results 
in §6. 

Throughout the paper we adopt the currently favoured 
cosmology Ho = 65 km s _1 Mpc -1 , Qm = 0.3, Q.a = 0.7. 



2 DLA OBSERVATIONS 

The data used in this paper are a compilation of several sur- 
veys at high (Storrie-Lombardi et al. f 996; Storrie-Lombardi 
& Wolfe 2000; Peroux et al. 2002) and low redshifts (Rao 
& Turnshek 2000; Churchill 2001). From these observations 
we extract information on three properties of DLAs which 
our models attempt to reproduce. 

The first, and most straightforward, quantity is the 
number density per unit redshift dn/dz — n(z), plotted in 
Figure 1. 

The second quantity is the column density distribution, 
obtained from the expression: 
n 

f(N, z)dNdX = =jj dNdX (1) 

where n is the number of DLAs with hydrogen column den- 
sity between N and N + AiV detected in the spectra of m 
QSOs encompassing a total absorption distance ^TLj AXi. 
The absorption distance is related to redshift by the expres- 
sion (Peebles 1993): 

x ( z ) = I Tr( 1 + *) 2 xE{z)dz (2) 
Jo H ° 

where 

e(z) = [n M (i + zf + n A }- 1/2 (3) 

Even with our relatively large sample, the statistics of 
f(N, z) dNdX are still limited, and necessitate dividing the 
data in only a few bins, both in AN and Az (see top pan- 
els of Figures 2 and 3). Also shown in these panels are the 
fits to the column density distributions derived by Peroux 
et al. (2002), assuming that f(N,z)dNdX is well approx- 
imated by a r function (similar to the galaxy luminosity 
function). These fits were based on a larger range of col- 
umn densities than that shown here, extending to the lower 
values of N(H I) appropriate to Lyman Limit systems (i.e. 
N(H I)< 2 x 10 2() cm" 2 ). 

The integral of the column density distribution gives us 
the third quantity of interest, the total mass of neutral gas 
traced by QSO absorbers. It is customary to express this as 
a fraction of the closure density: 

n DLA (*) = ^ / Nf(N,z)dN (4) 

Pcrit J Nm . n 

where thh is the mass of the hydrogen atom, /i is the mean 
atomic weight per baryon (fi = 1.4 for solar abundances; 
Anders & Grevesse 1989) and p cr it is the closure density 
(pcrit = 3#o/(8ttG) = 1.96 x 10~ 2 V gm cm" 3 , where h is 
the Hubble constant in units of 100 km s _1 Mpc -1 ). Values 
of Odla are plotted in Figure 4. 
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Figure 1. Redshift evolution of the number density of DLAs. 
The filled circles are data from Peroux et al. (2002) and refer- 
ences therein. The low redshift open circles are from Rao & Turn- 
shek (2000). The triangle at z = is from Churchill (2001). The 
lines show the number of absorbers predicted from the models 
described in §3. The solid line labelled 'TOTAL' is the prediction 
for all disk galaxies. The dotted and dashed lines show the respec- 
tive contributions to this total by high and low surface brightness 
disks. 

It is important to bear in mind that all the DLA sur- 
veys considered here used QSOs selected from magnitude 
limited samples. It has long been a concern (e.g. Pei & Fall 
1995) that in these surveys QSOs that lie behind dusty DLA 
systems will be under-represented, and that the statistics 
discussed above may therefore be incomplete and biased by 
an unknown amount. One way to determine the extent of 
such dust-induced bias is to consider the properties of DLAs 
in a complete, radio-selected, QSO sample, where dust is 
not an issue. The only such study conducted to date is 
the CORALS survey by Ellison et al. (2001). Although the 
CORALS statistics are still limited, it would appear that — 
at least as far as Odla is concerned — dust bias is a relatively 
minor effect, since the CORALS value is consistent with the 
earlier estimates within the errors of measurement (see lower 
panel of Figure 4) . 



3 MODELS OF DISK GALAXIES 

3.1 Basic Properties of the Models 

In Boissier & Prantzos (2000) a 'bivariate' family of spiral 
galaxies was described by considering disks with various ro- 
tational velocities, Vc, and spin parameters, A. In the local 
universe, the volume density Fv of galaxies with a given 
value of Vc can be deduced from the Tully-Fisher relation- 
ship and the luminosity function (Gonzalez et al. 2000). In 
the following, we adopt the distribution corresponding to 
the parameters listed in the fifth row of Table 4 of Gonzalez 
et al., with Vc between 80 and 360 km s" 1 . The spin param- 
eter A is a dimensionless measure of the angular momentum 
J of the dark matter halo surrounding the disk, and is given 
by: 

A = J^^G^M- 5 ' 2 (5) 
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Figure 2. Column density distributions at 2 < z < 2.7 (left), 
2.7 < z < 3.5 (centre) and 3.5 < z < 5 (right). In the top 
panels, the points are the binned observations and the dashed 
lines shows the best fits obtained by Peroux et al. (2002) for a T 
function extending to lower values of N(R I) than those shown 
here. The solid (dotted) lines show the model predictions at the 
lower (upper) boundaries of each redshift interval (i.e. at z = 
2.0 and 2.7 for the leftmost panel and so on). The lower panels 
show the differences between the model predictions and the T 
function fits to the distributions; the shaded areas show the ranges 
of uncertainty arising from the errors in the data (as shown in 
the top panels) and the ±0.3 dex typical uncertainty of the model 
predictions. 



where E and M are respectively the energy and mass of the 
dark matter halo. 

In the general framework of Cold Dark Matter theories, 
the final properties of the disk (its mass Md, scale-length 
Rd, and central surface density Eo) can be derived under 
simple assumptions on the relationship between the disk and 
the dark matter halo (see, for example, Mo, Mao, &White 
1998). The scaling relations 



Figure 3. Column density distribution at z < 2. As in Figure 
2, the top panel shows the binned observations of Rao and Turn- 
shek (2000) and the models for the redshifts z=0 and z=2. The 
dashed line is the power-law fit to the observations from Rao and 
Turnshck (2000). The lower panel shows the differences between 
the model predictions at redshift and 2 and the fit to the data, 
over the typical uncertainty as in Figure 2. 



(where the suffix MW refers to the Milky Way) are then 
used to relate these quantities to those of our Galaxy, as- 
sumed to be a typical spiral disk. 

Our models used these scaling relations to compute the 
chemical evolution in zones at different distances from the 
centres of disks with various rotational velocities and spin 
parameters. The history of each galaxy (for a given set of 
values of Vc and A) was derived by a 'backward' approach 
from observations of samples of nearby galaxies and cali- 
brated on the Milky Way (see Boissier & Prantzos 1999, 
2000; and Boissier et al. 2001 for more details). The resulting 
histories are characterised by an early formation of massive 
galaxies, while star formation occurs on longer timescales in 
less massive galaxies, following the infall of pristine gas. 
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3.2 The Spin Parameter Distribution and Low 
Surface Brightness Galaxies 

In our treatment we adopt the distribution of spin parame- 
ters derived from N-body simulations (Mo et al. 1998) and 
reproduced in Figure 5. Boissier & Prantzos (2000) consid- 
ered values for the spin parameter between A = 0.01 and 
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Redshift 

Figure 4. Redshift evolution of Hdla- The symbols have the 
same meaning as in Figure 1. As in Figure 1, the lines in the 
top panel show our model predictions for all disk galaxies and 
separately for high and low surface brightness components. In 
the lower panel, the solid line is the same as in the top panel, 
while the dashed line show the predictions of models without the 
'dust-filter' (see text). The filled triangle is from the CORALS 
survey by Ellison et al. (2001) which is free from the potential 
bias introduced by dust. 



0.09 and found that, within the framework of simple mod- 
els of chemical and spectrophotometric evolution of spirals, 
values of A > 0.07 apply to low surface brightness galaxies, 
with blue central surface brightness greater than 22.5 mag 
arcsec -2 , while galaxies with A < 0.07 are 'normal' spirals. 
As can be seen from Figure 5, the former account for 25% 
of the total number of galaxies in the distribution (integrat- 
ing between 0.07 < A < 0.21), while the latter contribute 
74% (integrating between 0.01 < A < 0.07). In this work we 
assumed that this population of LSB disks does indeed ex- 
ist and we accordingly extended the models of Boissier et al. 
(2001) to larger values of A than those considered in that pa- 
per in order to take fully into account this LSB population. 
Furthermore, we renormalised the distribution in Figure 5 
so that its integral over the range of values of A appropriate 
to normal, high surface brightness, disks (0.01 < A < 0.07) 
is equal to 1. This is because the normalisation of the dis- 
tribution of velocities Fv, derived from the observed galaxy 
luminosity function, does not include low surface brightness 
galaxies. While this treatment of LSB galaxies is simplistic, 
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Figure 5. Distribution of the spin parameter A. According to 
Boissier k. Prantzos (2000), galaxies with 0.01 < A < 0.07 are 
'normal' spirals, while galaxies with larger values of A are low 
surface brightness galaxies. The respective contributions of each 
(within the limits shown in the figure) to the total integral of the 
distribution are indicated. 



it does nevertheless allow us to test, in the framework of our 
model, their contribution to the DLA absorber population. 



3.3 Model Predictions for the Number Density, 
Column Density Distribution, and Total Gas 
Mass of DLAs 

Prantzos & Boissier (2000) and Hou et al. (2001) consid- 
ered the chemical evolution of the family of models of disk 
galaxies described above. In the zones within the galaxies 
where the gas density is sufficiently high to give a column 
density of neutral gas N(H I) > 2 x 10 20 cm -2 , the models 
match the chemical properties of DLAs, both in the overall 
degree of metal enrichment and in the more detailed abun- 
dance pattern of different chemical elements. However, in 
those treatments good agreement could only be achieved by 
imposing a 'dust filter', whereby DLAs with large column 
densities of metals (and therefore presumably dust) are ex- 
cluded from the samples, for the reasons explained in §2. 
This dust filter was based on the analysis by Boisse et al. 
(1998) and takes the form 



log (AT (HI) + [Zn/H]) < 21 



(9) 



where [Zn/H] is the abundance of Zn relative to solar (in the 
usual notation) and A r (HI) is the column density of hydrogen 
atoms (in units cm" 2 ). The filter operates in such a way that 
zones within a model galaxy which do not satisfy eq. (9) were 
not taken into account, on the assumption that such DLAs 
have been overlooked in existing surveys. 

Here we extend the work of Prantzos & Boissier (2000) 
and Hou et al. (2001) by using the models to make predic- 
tions for other properties of DLAs (specifically those dis- 
cussed in §2), with the aim of assessing what fraction of 
DLAs are the precursors of today's disk galaxies. We adopt 
the same 'backward' approach as in the earlier work and, 
like those authors, consider the effect of a dust filter defined 
as in eq. (9). Our treatment is broadly similar to that first 
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Figure 6. Surface and contour plots of the galactic cross-section (independent of how the mass is distributed between stars and interstellar 
gas) as a function of spin parameter A and rotational velocity Vc ■ 



adopted by Lanzetta (1993), but with the following major 
improvements: (i) we use models which predict the evolu- 
tion of present day galaxies back in time; (ii) we include 
in the models of a variety of galaxy morphologies, masses, 
sizes, and surface brighteness; and (iii) we are able to com- 
pare with the much larger data base of measurements of 
DLAs available now, ten years since the pioneering work of 
Lanzetta (1993). 

The number density of DLAs per unit redshift can be 
expressed as (Peebles 1993) 

»(*) = — = (- x -) = -(l + .) 2 x^)x G (.)(10) 

with E(z) as defined in eq.(3).* The quantity G(z) is the 
product of the volume density of galaxies and their cross- 
sections, averaged over the whole population 

G(z) = [ [ d 2 G (KV) (11) 
J \ Jv c 

where d 2 G(\,v) is the contribution to the total cross-section 
by galaxies with spin parameter A and rotational velocity 
Vc, given by 

d 2 G iKVc) = Fv(Vc)dVF x (X)dX / $dS-di (12) 

Js Ji=0 n 

Here J g dS represents the integral over the apparent surface 
of the galaxy, and J. di is the integration over all inclina- 
tion angles. The parameter $ is a switch which can take 
the values of 1 or 0. For example, it can be set to 1 inde- 
pendently of all other parameters in order to compute the 
geometrical cross section. The value of $ (1 or 0) can also 
be determined by the inclination and the local properties 
of galaxies (metallicity and gas density, both functions of 

* Note that the predicted value of therefore depends on the 
cosmological model assumed, whereas the observed value does 
not; the opposite is true for the column density distribution and 
^DLA- 



the distance from the galactic center) to exclude zones in 
galaxies which do not meet the DLA selection criteria, be- 
cause either of the conditions iV(H I)> 2 x 10 20 cm~ 2 or 
log (iV(HI) + [Zn/H]) < 21 (the dust filter in eq. 9) is not 
satisfied. Then, $ is set to in these zones, and 1 elsewhere. 

Figure 6 shows the function d 2 G(x,v c ) for $ = 1. From 
this figure it can be seen that, in our models, the major con- 
tribution to the cross-section for absorption (taken at this 
stage to be five times the scale length of the radial distribu- 
tion of mass without considering separately the distributions 
of the stellar and gaseous components) is from galaxies with 
relatively high rotational velocities (Vc — 200 km s _1 ) and 
intermediate values of spin parameter (A ~ 0.07), at the 
border line between high and low surface brightness galax- 
ies. This is straightforward to understand, since it is the 
product of Vc and A which determines the linear dimen- 
sions of galaxies. From the figure it would also appear that 
the ranges of values of Vc and A adopted are sufficiently 
broad to account for most of the cross-section. However, we 
may have underestimated the total cross-section if the faint 
end slope of the luminosity function is sufficiently steep that 
dwarf galaxies with Vc < 80 km s _1 make a non-negligible 
contribution. Such galaxies are not included in the models 
considered here because we are chiefly concerned with disk 
galaxies. 

We predict the column density distribution numerically 
by computing for each interval of column density [Nh, Nh + 
8Nh] the number of systems given by eqs. (11) and (12). 
The total neutral gas mass in DLAs is calculated from the 
expression 

n G = — — / F v {Vc)dVF x {X)d\ / $>Y, G dS-di, (13) 
Pc"t J\,v Js,i 71 

which is similar to eqs. (11) and (12), but with the second 
integral now inluding Eg, the mass surface density along the 
line of sight. The result of eq. (13) can then be compared 
with the observed value of Qola, as given by eq. (4). In 
both cases, $ was set to if N(R I) < 2 x 10 20 cm" 2 or if 
the condition imposed by eq. (9) was not met. 
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3.4 Neutral Atomic Gas 

The models of chemical evolution which we use to describe 
disk galaxies compute the surface density of the total amount 
of interstellar matter, in all its forms. On the other hand, 
for the comparison with the observed properties of DLAs we 
are interested only in the neutral atomic gas component. We 
do not consider the fraction of gas in ionised form because it 
is likely to be unimportant for our purposes. In the Galactic 
disk, for example, ionised gas accounts for less than 1% of 
the total gas mass of 3 x 10 9 M Q (Osterbrock 1974). Sim- 
ilarly, we can ignore the mass fraction in solid form (dust 
grains), which for solar composition is less than 2% (An- 
ders & Grevesse 1989) and would be even lower at the low 
metallicities typical of DLAs (e.g. Pettini et al. 1997). 

The correction for molecular gas is also likely to be of 
minor importance for the following reasons. From straight- 
forward geometrical considerations the cross-section for ab- 
sorption is always dominated by the outer regions of the 
disks where, in local galaxies at least, the molecular frac- 
tion is low. In the Milky Way, for example, H2 makes a 
significant contribution to the total gas profile only at small 
radii (in a molecular ring at ~ 4kpc), while the outer disk 
is totally dominated by atomic gas (e.g. Dame 1993). CO 
observations of nearby spiral galaxies show that the molecu- 
lar gas density decreases more steeply with radius than the 
atomic gas, and becomes a minor component beyond a few 
kpc (e.g. Boissier et al. 2002). In LSB galaxies the molec- 
ular fraction is observed to be very small (e.g. Matthews 
& Gao 2001). In any case, it is now well established that 
in DLAs iV(H 2 ) < JV(H I) (Petitjean, Srianand & Ledoux 
2002). Thus, in our models we take N(R I) = AT(Htot) (but 
correct the total gas density for the helium fraction). When 
we tried to compute the molecular fraction, using a recipe 
which links H2 to the star formation and calibrating on the 
Milky Way, we obtained results which are indistinguishable 
from those obtained assuming that all the gas is in atomic 
form. 



4 THE CONTRIBUTION OF DISK GALAXIES 
TO THE POPULATION OF DLAS 

We are now ready to compare the output of our models 
with the statistical properties of DLAs summarised in §2. 
The results are illustrated in Figures 1-4 and are discussed 
individually below. 

4.1 Number Density 

As can be seen from Figure 1, our models reproduce sat- 
isfactorily (bearing in mind the large uncertainties in the 
data) the observed number density of DLAs per unit red- 
shift in the interval < 2 < 2. It is interesting to note 
that our simplistic treatment of LSB galaxies, which essen- 
tially extrapolates the properties of normal spirals to larger 
values of the spin parameter A (see §3.2), predicts approxi- 
mately equal contributions to n(z) from high and low sur- 
face brighteness galaxies over this redshift interval. As we 
shall see, the available imaging data summarised in §5.1 are 
consistent with this rough breakdown. 

Figure 1 also shows that our models underpredict n(z) 



at z £ 2 by factors of several; while the data show a con- 
tinuous increase with look-back time, the predicted number 
density flattens off and decreases beyond z = 2. In reality, 
the observed increase in n(z) with z is consistent with the 
no-evolution scenario where the product of the comoving 
number density and cross-section of DLA absorbers remains 
constant in time — the measured n(z) increases simply be- 
cause the scale factor of the universe decreases with z. Thus, 
the flattening of n(z) in the calculations reflects the delayed 
and protracted epoch of disk formation in the models. From 
this we conclude that either today's disk galaxies were as- 
sembled earlier than predicted by our models, or a different 
class of objects is responsible for the majority of damped 
systems observed at z > 2. 

4.2 Column Density Distribution 

The observed and predicted column density distributions 
are compared in Figures 2 and 3. At z = 2 — 3 the agree- 
ment is fair, considering the large uncertainties in the data. 
The models tend to predict more high iV(H I) DLAs than 
observed, but this could be due to small number statis- 
tics (although it is interesting that a similar discrepancy is 
found with N-body simulations — see Katz et al. 1996; Schaye 
2001b; and Cen et al. 2002 for relevant discussions). At z < 2 
the agreement is also acceptable (Figure 3). Again we find 
a discrepancy at the high column density end of the distri- 
bution, but this time the models tend to underpredict the 
number of DLAs with logiV(H I) > 21.3 [N(R I) in cm~ 2 ] 
compared with the observations by Rao & Turnshek (2000) . 
It would clearly be highly desirable to improve the statistics 
at this end of the distribution with more extensive surveys 
for DLAs at all redshits. 

Returning to Figure 2, it appears that the agreement be- 
tween the models and the data becomes progressively worse 
as we move to z £ 3; the steepening of f(N) with increas- 
ing z uncovered by Peroux et al. (2002) is not reproduced 
by the models. Taken at face value, this is another indica- 
tion of a change in the population of galaxies responsible 
for DLAs; while the precursors of today's disks can account 
for the column density distribution at low and intermediate 
redshifts, at z £ 3 an additional population, with typically 
lower values of N(H I), seems to take over. 

4.3 Mass Density 

From Figure 4 it can be seen that our disk models under- 
predict Qdla by about a factor of two at all redshifts, al- 
though at z 2 the uncertainties in the determination of 
this quantity are so large that they do encompass our model 
predictions. 

All of the above results were obtained with the dust 
filter described in §3.3 (eq. 9). The lower panel of Figure 
4 shows the effect of removing such a filter and including 
in the statistics all DLAs, irrespectively of their metallic- 
ity and column density. We find that although J7dla does 
increase, the effect is not large, amounting to an upwards 
correction by a factor of 2 at z — 2 — 3, in good agree- 
ment with the empirical finding by Ellison et al. (2001). In 
retrospect this is not surprising. The condition imposed by 
cq. (9), log(iV(H I)) + [Zn/H]) < 21, is not satisfied primarily 
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Figure 7. H I sizes (radii at which jV(H I) = 2 X 10 20 cm~ 2 ) 
for face-on galaxies as a function of rotational velocity Vc for 
three different values of the spin parameter, A = 0.07 (solid line, 
considered here as the transition value between HSB and LSB 
spirals), 0.2f (dotted) and O.Of . (dashed). 'Normal' spirals occupy 
the horizontally hatched areas of the plot, while LSB galaxies fall 
in the vertically hatched regions). 



in massive galaxies (statistically less favoured by the shape 
of the luminosity function) and in the inner regions of galax- 
ies (statistically less favoured because they present a smaller 
cross-section than the outer regions). It is thus quite natu- 
ral that the dust filter should not remove a major fraction 
of potential DLAs. Figure 4 shows the dust bias being most 
important at z < 2. It remains to be seen, once the statis- 
tics of DLAs at these redshifts are improved, whether this is 
indeed the case. At present the uncertainties in the values of 
f^DLA deduced by Rao & Turnshek (2000) are so large that 
either set of models (with and without dust obscuration) is 
consistent with the observations. 



greater than the typical dimensions of galaxies today, both 
in their stellar Holmberg radii and in their cross-sections at 
N(H I) = 2 x 10 20 cm -2 . Figure 7 shows the cross-sections 
of galaxies in our models as a function of rotational velocity 
and spin parameter, at two different redshifts. Considering 
these plots in conjunction with those in Figure 6, it can be 
realised that such large radii are in fact not unexpected in 
our models once the contribution by LSB galaxies to the 
DLA population is taken into account. 



4.5 Comparison with Imaging Data at Low and 
Intermediate Redshifts 

In the last few years there have been many attempts to iden- 
tify the galaxies responsible for producing damped Lya sys- 
tems by deep imaging. We have collected in Table 1 all avail- 
able data from the literature pertaining to DLAs at z^ s < 1 
(where the results are most straightforward to interpret), re- 
stricting ourselves to cases where only one galaxy was found 
close to the QSO sight-line. Column (3) gives the projected 
QSO-galaxy separations in our Ho = 65 km s _1 Mpc , 

= 0.3, S1a = 0.7 cosmology. Column (4) and (5) list 
the absolute magnitudes of the putative absorbers in this 
cosmology, in either the rest frame B or K band; in de- 
riving these values from the observed magnitudes we ap- 
plied i\-corrections appropriate to the SED of an Sbc galaxy 
(Bruzual & Chariot 1993). The galaxy luminosities are given 
in units of L* in column (6), adopting AfjJ — 5 log h — —23.44 
(Cole et al. 2001) and M£ - 5 log h = -19.46 (Norberg et 
al. 2002). Column (7) gives an indication of the galaxy mor- 
phology as reported in the original sources referenced in the 
last column of the Table; cases where the identity of the ab- 
sorber has been confirmed with follow-up spectroscopy are 
indicated in column (8). 

The data collected in Table 1 form a small and far 
from homogeneous sample. Furthermore, in most cases there 
is no spectroscopic confirmation that the closest galaxy to 
the QSO sight-line is indeed the DLA absorber. Neverthe- 
less, these observations are all we have at our disposal at 
present, and it is still worthwhile to consider qualitatively 
how they compare with the predictions of our models. Con- 
sidering first the galaxy morphologies, the sample includes 
six apparently 'normal' spirals and five LSBs. This break- 
down is broadly in line with our expectations from Fig- 
ure 1 at z < l.T A further three galaxies do not appear 
to be spiral disks but have been classified by the original 
observers as irregular/dwarf /compact. The fraction may be 
higher because our sample in Table 1 does not include non- 
detections; nevertheless, preliminary indications are that 
this additional population of galaxies — which we have not 
considered here — is unlikely to make a major contribution 



4.4 Comparison with the Earlier Analysis by 
Lanzetta 

As mentioned above, our treatment is similar to that first 
originally implemented by Lanzetta (1993), although we 
have introduced several important refinements. It is thus 
of interest to compare our conclusions with his. In particu- 
lar, Lanzetta was forced to conclude that the incidence of 
DLAs required gaseous disks with typical radii of ~ 30 kpc, 



t It should be noted here that the definition of LSB in our mod- 
els is based on whether the spin parameter A is larger than 0.07 . 
Owing to the larger star formation rate in the past in the models, 
some galaxies with A > 0.07 at z ~ 1 actually have values of cen- 
tral surface brightness somewhat greater than 22.5 mag arcsec -2 . 
However, this is not expected to affect significantly the finding 
that spirals and LSB galaxies make comparable contributions 
to the DLA cross-section because the morphology classifications 
listed in Table f are only approximate estimates. 
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QSO 


2 


D (kpc) 


M B 


M K 


L/L, 


Type 


Confirmed 


Refs 


HS1543+5921 


0.009 


0.5 


-15.7 


-19.3 


0.01 


LSB 


yes 


(a) 


OI363 


0.0912 


< 3.7 




-22.1 


0.12 


LSB 




(b) 


OI363 


0.2213 


23 




-22.4 


0.16 


Dwarf 


yes 


(b) 


0952+179 


0.239 


8.1 




< -20.1 


< 0.02 


Patchy/irr LSB 




(b,g) 


1127-145 


0.313 


19.7 




-20.9 


0.04 


Patchy /irr LSB 


yes 


(b,g) 


PKS1229-021 


0.3949 


8 


-19.8 




0.6 


Irr. 




(c,d) 


3C196 


0.437 


9 


-21.5 




2.8 


Spiral 




(c) 


0827+243 


0.518 


40.2 


—21.3 


—25.3 


2.3 


Spiral? 


VPS 


(b s) 


Q1629+120 


0.532 


20.3 


-20.8 


-24.8 


1.5 


Spiral 




(g) 


0058+019 


0.6118 


8.7 


-19.8 




0.6 


Spiral 




(e) 


Q1209+107 


0.6295 


11.9 


-21.1 




1.9 


Spiral 




(c) 


3C286 


0.692 


19.1 


-20.7 




1.3 


LSB 




(d) 


MC1331+170 


0.7443 


30.4 


-21.6 




3.0 


Edge-on Spiral 




(c) 


PKS0454+039 


0.8596 


6.6 


-19.9 




0.6 


Compact 




(c) 



Table 1. DLAs and candidate galaxies. References: (a) Bowen et al. (2001); (b) Turnshek et al. (2001); (c) Le Brun et al. (1997); (d) 
Steidel et al. (1994); (e) Pettini et al. (2000); (f) Steidel et al. (1997); (g) Nestor et al. (2001). 




20 40 60 

Impact Parameter (kpc) 



Figure 8. Distributions of impact parameters. Top panel: distri- 
bution constructed from Table 1 grouping the data in 10 kpc bins. 
The white areas of the histogram correspond to galaxies classified 
as irregular /dwarf /compact. Middle and bottom panels: distribu- 
tions predicted by our disk models at z = 1 and 0. 



to the DLA statistics. Only a uniform survey to well denned 
detection limits can address this question quantitatively. 

Figure 8 compares the observed and predicted distribu- 
tions of projected impact parameters D; the latter were cal- 
culated for random orientations of the galaxies on the plane 



of the sky. Within the limitations imposed by the small size 
of the sample, the observations apparently favour smaller 
values of D than predicted. In some cases, misidentifications 
of the closest galaxy to the QSO sight-line as the absorber 
may be the reason for this apparent discrepancy. The models 
predict an extended tail at large impact parameters which 
cannot yet be discerned with the limited data available; it 
will be of interest to test this prediction with future surveys. 

Figure 9 compares the observed and predicted distribu- 
tions of Mb ■ To improve the limited statistics, we have used 
the values of L/L* in column (6) of Table 1 (which in all but 
three cases were based on either B-band or if-band magni- 
tudes, but did agree in the three cases where magnitudes in 
both bands are available) to calculate M b from the relation 
Mb - 5 log ft = -19.46 (Norberg et al. 2002). Evidently, our 
models of disk galaxies can reproduce reasonably well the 
magnitudes of the brighter galaxies associated with DLAs, 
but do not extend to the low luminosities of the faintest ab- 
sorbers in the sample. Again, the relatively high frequency of 
faint LBGs in the present sample of DLA galaxies needs to 
be investigated further with more extensive imaging surveys 
before firm conclusions can be drawn. 



5 THE POPULATION OF DLAS AT HIGH 
REDSHIFTS 

As discussed in §4 and shown in Figures 1-4, our models for 
disk galaxies provide a progressively poorer match to ob- 
servations of DLAs beyond z ~ 2 — 3. Taken at face value, 
this may indicate that the precursors of present day spi- 
ral galaxies — evolved back in time according to our models 
which do not include merging — are insufficient to account for 
the observed numbers of DLAs and that an additional pop- 
ulation of absorbers is required. We investigate this point 
further in Figure 10. In the top panel we show, on a log 
scale and as a function of redshift, the ratio of the observed 
number of DLAs (per unit redshift) to the number which 
can be attributed to the progenitors of today's disks accord- 
ing to our models. Thus this plot shows the ratios between 
the data points and the line labelled 'TOTAL' in Figure 1. 
While at z ~ 2 disks can account for the number density 
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Figure 9. Distributions of absolute B-band magnitudes. Top 
panel: distribution constructed from Table 1 grouping the data in 
1 magnitude bins. The white areas of the histogram correspond 
to galaxies classified as irregular/dwarf/compact. Middle and bot- 
tom panels: distributions predicted by our disk models at z = 1 



and 0. The vertical dotted line at 
of Ait, in our cosmology. 



-20.4 indicates the value 



of DLAs within the uncertainties in both the data and the 
model predictions, at redshift z ~ 4.5 DLAs are nearly ten 
times more numerous than predicted. This excess evidently 
decreases rapidly with time, approximately by a factor of 
two per unit redshift in the interval 2 < z < 4.5 (dashed 
line) . 

The bottom panel of Figure 10 shows the contribution of 
our model disk galaxies to the observed value of Hdla (the 
ratios between the data points and the line labelled 'TO- 
TAL' in the top panel of Figure 4). Also shown in this panel 
is the progressive transfer of baryons from the overall gas 
reservoir — from which disks are formed in our models — into 
galaxies, as disks form and grow with time. Evidently, disks 
can account for approximately half of J7dla at all epochs 
(as already shown in Figure 4). In other words, it appears 
that the additional population of absorbers which dominate 
the number counts of DLAs at the highest redshifts do not 
make a significant contribution to £7dla compared with the 
baryons which are in disks. This is reflected by the steepen- 
ing of the column density distribution with increasing red- 
shift found by Peroux et al. (2002) and shown in Figure 2. 

A possible interpretation of these results is that the 'ex- 



Figure 10. Top: The ratio between the observed number den- 
sity of DLAs per unit redshift and the value predicted by our 
disk galaxy models. The data suggest the existence of a rapidly 
evolving 'excess' population of DLAs at high z over what can be 
accounted for by extrapolating the properties of today's spiral 
galaxies back in time. The dashed line corresponds to a halving 
of n(z) for this additional population per unit redshift between 
z = 4.5 and 2. Bottom: Ratio between the mass density of the 
total gas reservoir available for disk formation and that which has 
been incorporated into disks (dashed line); and between the ob- 
served mass density in DLAs and that in our model disks (filled 
circles). 



cess' population of DLAs at high redshift are small struc- 
tures, of relatively low column density, which progressively 
merge with each other and with more massive galaxies to 
form the progenitors of today's disk galaxies by z ~ 2. 
This scenario, which at the moment can only be regarded 
as speculative, would explain at the same time the observed 
decrease of n(z) and flattening of the column density dis- 
tribution (resulting in relatively more numerous systems of 
high column density) with the progress of time from z = 4.5 
to 2. Observationally, this scenario is consistent with: (a) 
the initial results of imaging searches for the galaxies re- 
sponsible for DLAs at 2 ~ 4 which suggest that they are 
fainter than ~ 1/4L* (Prochaska et al. 2002); and (b) the 
low metallicity — typically 1/100 of solar, only marginally 
higher than that of the intergalactic medium traced by the 
Lya forest — of DLAs at these high redshifts (Prochaska & 
Wolfe 2002). 
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6 SUMMARY AND CONCLUSIONS 

In this paper we have used the models developed by Boissier 
& Prantzos (2000), which describe the chemical and spec- 
trophotometric evolution of disk galaxies of different surface 
brightnesses, to assess whether the progenitors of present 
day spirals can account for the population of DLA absorbers, 
as is often assumed. To this end, we have compared the 
model predictions with several, recently determined, statis- 
tical properties of DLAs, in particular their number den- 
sity per unit redshift, their column density distribution, and 
its integral which gives Odla, the total mass of neutral gas 
traced by DLAs. In addition we have brought together avail- 
able imaging data on these absorbers at z < 1 from the 
literature, to compare with the morphologies, sizes and lu- 
minosities of the model galaxies. Our main results can be 
summarised as follows. 

1. The models are reasonably successful at reproducing 
many of the properties of DLAs at redshifts z 2, within the 
uncertainties in the measurements. Specifically, the number 
density of absorbers per unit redshift and about half of the 
neutral gas mass they trace can be reproduced with evolving 
disk galaxies in our models. Dwarf galaxies, which have not 
been considered here, may account for some of the DLA 
absorbers, but this population is not the dominant one. 

2. Normal spirals and low surface brightness galaxies 
make comparable contributions in our models to both the 
numbers of DLAs and their neutral gas mass. Turning this 
statement around, it is only with the inclusion of LSBs that 
it is possible to reproduce the statistics of DLAs in our 
models. While LSBs may not be as numerous as high sur- 
face brightness galaxies, their large dimensions and high gas 
content combine to make a significant contribution to the 
overall cross section for H I absorption. The scant imag- 
ing data available at z < I is broadly in agreement with 
this conclusion, although the observed distribution of im- 
pact parameters D seems to be narrower, and more peaked 
towards lower values of D, than our models predict. The 
distributions of luminosities agree at the bright end, but the 
data also include a few cases of very faint DLA-producing 
galaxies which are not predicted by the models. More uni- 
form and extensive imaging surveys are required to reach 
firm conclusions on both of these points. Furthermore, the 
properties of LSBs, which have been included in our models 
according to a rather simple recipe, need to be investigated 
further for a proper comparison with DLAs. 

3. We have investigated the effects of a possible dust- 
induced bias in current DLA samples, which may lead to an 
underestimate of the relative number of sightlines through 
galactic regions where the column densities of gas and metals 
are both high. We have found this potential problem to have 
a relatively minor effect, particularly on Odla, in agreement 
with the initial results from the CORALS survey by Ellison 
et al. (2001). 

4. As we move from z ~ 2 to higher redshifts, models 
based on extrapolating back in time the properties of today's 
disk galaxies are progressively less successful in reproducing 
the statistics of DLAs. The data can be interpreted as evi- 
dence for the existence of an additional population of DLAs, 
of generally lower column density, which dominate the num- 
ber density of absorbers over the progenitors of today's disks 
beyond z ~ 3. Possibly these are sub-units destined to merge 



with each other and eventually with more massive galaxies 
by z ~ 2. This interpretation is supported by the failure 
up to now to detect the galaxies producing DLAs at z ~ 4 
in deep images reaching down to 1/4L*. The history of 
assembly of today's galactic disks may well be reflected in 
the redshift evolution of the number density of DLAs. 

Samuel Boissier would like to acknowledge the support 
of a Framework 5 Marie Curie fellowship through contract 
number HPMF-CT-2000-00521. 
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